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ABSTRACT 

The p Cephei stars represent an important class of massive star pulsators probing the evolution of B-type stars and the transition 
from main sequence to hydrogen-shell burning evolution. By understanding p Cep stars, we gain insights into the detailed physics 
of massive star evolution such as rotational mixing, convective core overshooting, magnetic fields and stellar winds, all of which 
play important roles. Similarly, modeling their pulsation provides additional information into their interior structures. Furthermore, 
measurements of the rate of change of pulsation period offer a direct measure of p Cephei stellar evolution. In this work, we compute 
state-of-the-art stellar evolution models assuming different amounts of initial rotation and convective core overshoot and measure 
theoretical rates of period change for which we compare to rates previously measured for a sample of p Cephei stars. The results of 
this comparison are mixed. For three stars, the rates are too small to infer any information from stellar evolution models, whereas for 
three other stars the rates are too large. We infer stellar parameters, such as mass and age, for two p Cephei stars: CMa and S Cet, 
that agree well with independent measurements. We explore ideas for why models may not predict the larger rates of period change. 
In particular, period drifts in p Cep stars can artificially lead to overestimated rates of secular period change. 
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1. Introduction 

Massive B-type stars are powerful engines which drive cosmic 
evolution from star formation to galaxy evolution to chemical 
enrichment and the seeds of life in the Universe. These stars im¬ 
pact star formation through their winds contributing momentum 
and turbulence into the interstellar medium and chemically en¬ 
rich the ISM through their supernovae. Exploring the physics 
of these massive stars contributes to the understanding of cos¬ 
mic evolution as well as stellar physics in general. Massive stars 
are also important probes of the physics of stellar evolution be¬ 
cause these B-type stars display many different kinds of phe¬ 
nomena related to the physics or stellar rotation, stellar winds, 
pulsation, binarity, magnetic fields, and possibly stellar merg¬ 
ers ( Langer| 2012| l. Eor instance, the emission line stars dis¬ 
play rapid, nearly critical, rotation, while the fi Cephei stars 
undergo radial and non-radial pulsation ( Townsend et al.|2004[ 
|Moskalik & Dziembowski|1992| l. In this work, we will explore 
the p Cephei variable stars because their pulsation allows us to 
study their internal structure and potentially their evolution. 

The p Cephei variable stars are particularly interesting be¬ 
cause they display a diverse range of phenomena. Some p Cep 


stars appear to have strong magnetic fields (Silvester et al.|2009| 

|Hubrig et al.||2011| |Neiner et al.||2012|l while others have weak 

or no magnetic fields (|Silvester et al.|2009 

lEossati et al. 

2015 1, 

some rapidly rotate (lAerts et al. 1994 E 

andler et al. 

20121, 

while others rotate slowly (IShultz et al. 

20151. Most, if not 


2005 [ Huenemoerder et al. 20121 and are also X-ray sources. 
Recently, Oskinova et al. ^ 2014|l discovered X-ray variability in 


one p Cephei star CMa over a pulsation cycle. Understanding 
the physics in these stars impacts our understanding of massive 
star physics in general along with the evolution from the main 
sequence to the red supergiant stage of evolution and beyond. 


While probing the physics of massive stars, the p Cephei 
stars also span the boundary between the main sequence evo¬ 
lution and the transition to the blue supergiant stage of evolution 
( Demarque & Percy] 1965 [ Stothers|1965 1. As such, observations 
and models of these stars help us understand the transition and 
precisely where on the Hertzsprung-Russell diagram it occurs 
thus, again, providing opportunities to refine the input physics in 
stellar evolution models. 

In particular, there are two aspects of stellar physics that 
can be tested using models of p Cephei variable stars and ob¬ 
servations: convective core overshooting and rotational mixing 
( |Miglio et al.||2009) [Lovekin & GoupiT]|20101 l. Convective core 
overshooting mixes material into the stellar core from layers 
some fraction of a pressure scale height above. In stellar evo¬ 
lution models, mixing length theory treats convection as being 
fixed within the convective region such that the velocity and ac¬ 
celeration of a convective eddy goes to zero at the boundaries. 
Convective core overshooting is an ad hoc prescription in stellar 
evolution codes designed to account for the fact that convective 
eddies will not have zero velocity at the boundary between con¬ 
vective and radiative regions, hence penetrating above the con¬ 
vective core and mixing additional hydrogen into the stellar core. 
This overshooting acts to extend the main sequence life time of a 
massive star by create a more massive core at the end of main se¬ 
quence evolution. Rotational mixing will have similar effects on 
massive star evolution by meridional mixing (e.g. |Erischknecht 
let al.|2010trBrott et aLpOlTbl l. 

Thanks to asteroseismic observations of p Cephei variables, 
there have been numerous measurements of the amount of con¬ 
vective core overshoot, parameterized as a fraction of a pres¬ 
sure scale height, acoHp, as the total amount of material mixed 
over a main sequence lifetime. The typical amount of convective 
core overshoot required to match observations is about a^o = 0.2 
and is consistent with measurements for the p Cep stars p CMa 
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( [Mazumdar et al.||2006| l and 6 Ceti (|Aerts et al.l|2006[) as well 
for classical Cepheids ([Cassisi & Salaris|201 if Neilson et al. 


2011 Neilson & Langer [2012 Prada Moroni et al.||2012 i and 

other stars. However, measurements for other fi Cep stars require 
different amounts of overshoot such as 6 Oph a^o = 0.4 (|Briquet 
et al |2007|l 0.28//p (|Lovekin & Goupil|20T0l l, HD 129929 


0.1 (|Aerts et al.|2003| l, and v Eridani Oco - 0.28 ( jSuarez et al.| 


2009|l. These measurements are obtained by comparing pulsa 


tion and evolutionary models of p Cep stars. 

The degeneracy between rotational mixing and convective 
core overshooting is difficult to break. Rotational mixing can 
be constrained from measurements of the [N/C] enhancement 
( [Hunter et al.|2008]|Brott et al.|201 lb} , but this enhancement has 
been detected in some slowly-rotating stars ( [Morel et al.|2006| l. 
There is no obvious way to constrain the importance of one pro¬ 
cess relative to another. The goal of this work is to compare new 
stellar evolution models with rotation to rates of period change 
measured for Galactic p Cephei stars and test whether period 
change is a probe of rotation and overshoot. 

Eggleton & Perc}^ ( 1973|l first compared stellar evolution 


models with period change measurements and determined that 
the rate of period change is sensitive to the stage of evolu¬ 
tion, be it main sequence or blue supergiant. However, since 
that work, stellar evolution models have matured, incorporat¬ 
ing new opacities plus improved prescriptions for physics of 
stellar mass loss, rotation and magnetic fields (Yoon & Danger 


2005 [ l. For example, these changes have led to the realization that 
P Cephei pulsation is driven by the iron opacity bump (Moskalik 
& Dziembowski[[l992[ [Cox et ar][1992[ [Dziembowski et ah 
1993\ as opposed to helium ionization in classical Cepheid and 


RR Lyrae stars. With the advances in stellar evolution models, 
it is an opportune time to revisit the comparison with period 
change measurements. 

In the next section, we discuss the stellar evolution models 
used in this work, and compute theoretical rates of pulsation pe¬ 
riod change in Sect. 3. In Sect. 4, we compare the rates of period 
change and test how period change relates to rotation rates. In 
Sect. 5, we discuss and summarize the results. 


2. Models 


We compute stellar evolution models using the Yoon & Danger 
( 2005| l code for masses ranging from M = 7 to 20 Mq in steps 
of one solar mass with various assumptions for internal mixing. 
One grid of models is computed assuming moderate convective 
core overshooting at Oco = 0.2 (as described by for example 
Huang & Weigert 1983[[Neilson et al. 20111 and with initial rota¬ 
tion rates Vrot = 100,200 and 400 km s V We include rotational 
mixing in these models based on the prescription of[Heger & 


Danger 


(2011a 


(2000 1 , which is described in more detail by Brott et al. 
A second grid of models is computed without rotation 


but assuming convective overshooting parameters = 0,0.2, 
and 0.4. 

The stellar evolution models are computed assuming 
radiatively-driven winds as prescribed by Vinket aL| (2000l dur¬ 
ing main sequence and blue supergiant evolution and with a 
Grevesse & Sauval ( 1998| l solar metallicity. We compute the 
evolution models from the zero-age main sequence past the 
boundaries of the p Cephei instability strip as measured by 
Pamyatnykh ( 2007[ ). Sample stellar evolution tracks are plotted 
in Fig. [IJ for different rotation rates and convective core over¬ 
shooting parameters. 


Table 1. Rates of period change for a sample of p Cep stars from 


Jerzykiewicz (1999 1 . 


Star 

Period (d) 

PjP (Myr-‘) 

7 Peg 

0.15175 

+0.046 ±0.015 or 0.0 ±0.11 

(5Cet 

0.16114 

+0.20 ± 0.01 or 0.34 ± 0.06 

CMa 

0.20958 

+0.20 ± 0.03 

a Lup 

0.25985 

+0.84 ± 0.22 

cr Sco 

0.24684 

+1.55 ±0.14 

V2052 Oph 

0.13989 

-0.10 ±0.22 

BW Vul 

0.20104 

+1.35 ±0.02 

P Cep 

0.19049 

+0.06 ± 0.06 


3. Period Change 

The purpose of this work is to compare measured rates of pe¬ 
riod change from observations with those predicted from stellar 
evolution models representing p Cephei stars. To do so, we use 
an analytic mod el of period chan ge based on the period-mean 
density relation ([Eddington| 191 8|l, P ^[p - Q, similar to that of 
Eggleton & Percy|([19 /3|l. We do not use the derivations from 


Neilson et al.[ P0i2a|b|l because there is no evidence that the 


pulsation constant, Q varies as a function of period for p Cephei 
variables. Therefore, we write the relative rate of period change 
as 


P 

P 


IM 3R 
'2M 2R’ 


( 1 ) 


where M the rate of change of stellar mass. This relation mea¬ 
sures the evolution of the star as a function of time. 

We computed relative rates of period change from the stel¬ 
lar evolution models using this relation. Rates of period change 
are shown in Fig. |^for two sets of models. The hrst set of mod¬ 
els are computed assuming zero initial rotation, but with differ¬ 
ent amounts convective core overshooting, while the second set 
of models assume moderate convective core overshooting only 
but varying the initial rotation rates. This result suggests that 
the combination of measured rates of period change and a fun¬ 
damental stellar parameter such as stellar radius, effective tem¬ 
perature or luminosity cannot provide a meaningful measure of 
physical processes such as convective core overshooting and/or 
rotational mixing (e.g. [Miglio et al.[|2009] l. These mixing pro¬ 
cesses act to change the mean structure of the star, but only in a 
way that is degenerate with models with different stellar masses 
but no additional mixing. 

As part of this analysis, we use rates of period change 
presented previously by Jerzykiewicz & Pigulski| ( \999\ and 
[Jerzykiewicz[ \ \999\ for eight p Cephei stars, see Table [T] It 
should be noted that these rates differ from those presented 
by Eggleton & Percy ( [1973| l, based on different analyses. 
Measuring period changes for p Cephei stars is a challenge as 
for many of these stars there are multiple oscillation frequencies 
and the periods might confuse the measurement. A second issue 
has been addressed by [Pigulski[ ( [1992[ [1993[ [1994[ l, for which 
binary companions affect measurements of period change due to 
the light time effect. While this is not a complete sample of pe¬ 
riod change measurements for p Cephei stars, it does represent a 
well-measured sample. 


4. Results 

We compare theoretical and measured rates of period changes 
for the sample of p Cep stars in Table [T] From these measure- 
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Fig. 1. (Left) Stellar evolution tracks assuming moderate convective core overshoot, = 0.2 and initial rotation rates of 0,100,200, 
and 400 km s“' denoted by black solid, green dashed, blue dot-dashed and magenta dotted lines, respectively. (Right) Stellar 
evolution tracks assuming zero rotation and convective core overshoot, Oc - 0,0.2, and 0.4 as denoted by black solid, green dashed 
and blue dot-dashed lines, respectively. 


ments, we can compute stellar fundamental parameters for mod¬ 
els with rates of period change consistent with each observation. 
Based on the rates presented in Tableand the models plotting 
in Fig.[^ we break the comparisons into three categories; 1) pe¬ 
riod change rates fit by almost all models, 2) period change rates 
fit by specific models and 3) period change rates fit by no mod¬ 
els. The first category simply suggests that the measured rate 
of period change offers no constraints on the models, while the 
third category suggests that the period change is completely in¬ 
consistent with models. These differences appear to be a result 
of insufficient precision for some observations and that obser¬ 
vations also span a greater range of period changes rates than 
predicted from the stellar evolution models. 


There are three (i Cep stars for which the measured rates of 
period change offer no constraints: y Peg, V2052 Oph, and the 
prototype (3 Cep. These three stars have both the smallest val¬ 
ues for the rate of period change but also with the greatest error 
bars. Similar, the stars that fall into the third category include 
those with the largest rates of period change: a Lup, cr Sco, and 
BW Vul. While the | Jerzykiewicz | (19991 sample consists of only 
eight stars (nine measured rates of period change), there are sig¬ 
nificant hints that there is something amiss in stellar evolution 
models that was not found by Eggleton & Percy] ( |1973[ ), even 
though their models did predict large rates of period change. 
This is because their models predict that stars evolve as blue 


supergiant stars within the yS Cephei instability strip. All of the 
P Cephei stars in the sample were classified as either luminosity 
class III or IV. The three p Cep stars with large rates of period 
change are luminosity class III but so is the prototype p Cep 
that has a measured rate of period change consistent with zero. 
There is no obvious connection between rates of period change 
and evolutionary state of the stars. 


Is there anything special about the three p Cephei stars: 
a Lup, cr Sco, and BW Vul? These three stars appear to pulsate 
consistently in the fundamental mode, but there is nothing appar¬ 
ently unique about these stars relative to the others in the sam¬ 
ple. The mystery deepens when one notes that the other p Cephei 
stars in our sample tend to have magnetic fields, which are not 
modelled in this work, and many are slow rotators. The stars that 
are consistent with none of the stellar evolution models are the 
ones that are expected to be most easily fit. 


The three stars with rates of period change that are consis¬ 
tent with predictions from all stellar evolution models, y Peg, 
V2052 Oph, and p Cep, also provide no clear answers. There 


have been claims that y Peg has a weak magnetic field 

(Butkovskaya & Plachinda 

2007|l, but Neiner et al. (2014)1 re- 

ported a null measurement. 

’he star also has a measured rotation 

rate of about 3 km s ' (jTelting et al. 2006 

Handler 2009 

1 . These 


observations suggest the stellar evolution models computed in 
this work should be sufficient for analyzing the stellar proper- 
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Logio(L/L0) 


Logio(L/L0) 


Fig. 2. (Left) Rates of period change computed from stellar evolution tracks assuming moderate convective core overshoot, = 0.2 
and initial rotation rates of 0,100,200, and 400 km s“^. (Right) Rates of period change computed from stellar evolution tracks 
assuming zero rotation and convective core overshoot, ac - 0,0.2, and 0.4. The different colors are same as for Fig. [T] while the 
horizontal orange line denotes the measured rates of period change for CMa and 6 Cet and the gold line denotes the second 
measured value for i5 Cet. 


Table 2. Measured fundamental properties of two p Cephei stars 
from the literature. 


Property CMa“ 6 Cet'’ 

Mass (Mo) 13-15 T9 

T,f f{K) 25900 ± 100 21675 ±400 

Radius (Rq) 8.7 ± 0.7 — 

logL/Le 4.46 ±0.07_3.6 


Note 1. (a) |Shultz et al.| ( [2015| l, (b) |Levenhagen et'aL] ( |2013) 


ties. On the other hand, V2052 Oph is a magnetic He-strong star 
(Neiner et al. 2012[ ). Briquet et al.| ( 2012[ ) measured rotation rates 
of about 70-75 km s ' and smaller amounts of convective core 
overshooting than for other p Cephei variable stars. The proto¬ 
type p Cep is similarly a magnetic He-strong star with a varying 
magnetic field of about 100 G ( Henrichs et al.|20T3 l. The three 
P Cep stars that cannot be constrained, i.e., consistent with all of 
the computed models, appear to have little in common. 

Of the two stars for which we find rates of period change that 
can measure stellar properties: CMa and 6 Cet, the former has 
a measured magnetic held strength of about 100 G (Hubrig et al. 
2006 Silvester et al.|2009)l and very slow rotation. This behavior 


is consistent with spin down due to the magnetic held interac¬ 
tions ( [Shultz et al.|2015) . There have been no measurements of 
magnetic helds in 6 Cet, but it does have a nitrogen excess and 
slow rotation rate consistent with known magnetic p Cep stars 
( Morel et al.|2006[ Levenhagen et al.|2013| l. 

In the next subsections, we explore in more detail the pre¬ 
dicted stellar properties for f' CMa and 6 Cet in combination 
with measurements of rotation rates and measurements of fun¬ 
damental parameters from seismic and other analyses. 


4.1. Properties of CMa 


Shultz et al.j ( [2015 1 presented a Baade-Wesselink analysis of 
CMa using Hipparcos photometry and radial velocity mea¬ 
surements (jSaesen et al.|2006|). They measure a mean effective 


temperature of 25.9 + 0.1 kK, a radius R = 8.7 ± 0.7 Rq and 
a luminosity log LILq = 4.46 + 0.07. The authors also mea¬ 
sure the age of CMa to be about 12.6 Myr from isochrones 
( Ekstrom et al.|20f2 i. Furthermore, they note that the star is ro¬ 
tating slowly with a period of about 60 years. 

Using the measured rate of period change and our stellar evo¬ 
lution models assuming zero rotation only we measure an age 


between 8-18 Myr, consistent with results from the Ekstrom 
jet al.| ( |20T^ stellar isochrones. However, when we include the 
effective temperature as an additional constraint of the stellar 
parameters, we measure a luminosity between logL/Lo = 4.40 
- 4.60, but with a younger age of 8 - 12 Myr. Eurthermore, we 
measure the mass of CMa tobeM = 13-15 Mq with a small 
probability of it being more massive if there is no convective 
core overshooting in the star. 


This comparison is striking because of how consistent the 
measured rate of period change is with predictions from our 
models while for other cases predicted rates are too small rel¬ 
ative to observations. Another issue is that CMa has a strong 
magnetic field, which is not included in our models. If the mag¬ 
netic field is strong enough then one might expect it to affect 
the pulsation properties of the star, but these findings suggest the 
magnetic field plays no role in the period change. It is not clear 
why the models provide such consistent predictions of stellar pa¬ 
rameters for this star but not for others ([Mathis & Neiner|2015]l. 


4.2. Properties of 6 Cet 


The second case where predicted rates of period change ap¬ 
pear to constrain stellar fundamental parameters is for b Cet. 
Levenhagen et al. (2013 1 fit spectroscopic observations with 
non-LTE models and measured T^ff - 21675 + 400 K and 
logg = 4.03 + 0.08. Other reports measure similar values for 
the effective temperature but gravities as small as logg = 3.7 


(jNiemczura & Daszyhska-Daszkiewicz 

2005l[Morel et al.[2006 

[Hubrig et al.ll2009|r 

Levenhagen et al. 

(2013|l also measured a 

projected rotation ve 

ocity of about v sin / = 13 km s ', suggest- 


ing that the star is slowly rotating or is being viewed nearly pole- 
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on. Comparing their results to stellar evolution models, they also 
infer a mass = 7.9 Mq and a luminosity of logL/Lo = 3.6. 

When we compare our stellar evolution models to both mea¬ 
surements of period change for 6 Cet, we find weakly con¬ 
strained values of stellar parameters. The predicted age, from 
7-18 Myr, is consistent with previous results ([Levenhagen 


et al.|2013]l. The effective temperature is predicted to range from 


logTeff = 4.32 - 4.42, and the luminosity is logL/Lo = 4.2 
5. The effective temperature is consistent with observations but 
the luminosity is significantly greater. There is no improvement 
if we consider only stellar evolution models with no rotation as 
opposed to the case of CMa. 

If we confine our models to only those with effective tem¬ 
peratures consistent with the [Levenhagen et'af ( 2013| l measure¬ 
ments, we better constrain the other stellar parameters. In this 
case, we measure the mass of 6 Cet tobeM= 10-12 Mq 
and its luminosity is logL/Lg = 4.18 - 4.25. The luminosity 


is still significantly greater than that measured by Levenhagen 
|et al.| ( |20T3| ) using [Schaller et ar] ( |1992[ ) stellar evolution tracks. 
The differences could be due to different prescriptions of interior 
mixing, opacities and input physics, making it unlikely that our 
results indicate some unknown underlying issue with the stellar 
evolution models. 

While a magnetic field has not been directly measured for 
6 Cet, the star shows nitrogen enhancement which is suggestive 
of one ( Morel et al.|2006) Levenhagen et al.|2013| l. When cou¬ 
pled with the apparent slow rotation of 5 Cet, it appears likely 
that 6 Ceti is similar to CMa. This result further raises ques¬ 
tions as to why stellar evolution models appear to fit these two 
stars, but not the non-magnetic fi Cephei variables in the sample. 


5. Discussion 


are points in the evolution of every star along the p Cephei in¬ 
stability strip where the rate of period change is very small. For 
three other stars: BW Vul, cr Sco, and a Lup, the predicted rates 
of period change for all stellar evolution models are too small 
relative to the measured rates. For two stars CMa and 6 Cet, 
we predicted fundamental stellar parameters by comparing pre¬ 
dicted and observed rates of period change. The predicted stellar 
parameters are consistent with previous independent measure¬ 
ments. We also measure the masses for d Cet and x/' CMa to be 
M - \ \ ±\ Mq and M - \A±\ Mq, respectively. 

The fits to the measured rates of period change of CMa 
and 6 Cet are both reassuring and surprising. That we find agree¬ 
ment with our models suggests that the analysis is consistent and 
the models are good representations of real stars. However, there 
is evidence for strong magnetic fields in the two stars, but these 
are not taken into account for the stellar evolution models. That 
the models tend to agree suggests that the presence of a strong 
magnetic field might not play a significant role in the pulsation 
of these stars relative to secular stellar evolution. 

Those stars with measured rates of period change that are 
greater than any of those predicted by stellar evolution models 
is disconcerting. No amount of rotation or convective core over¬ 
shooting can resolve that difference nor do these three stars have 
any common property that could explain the difference. For in¬ 
stance, only V2052 Oph shows evidence for a strong magnetic 
field. There are a number of reasons as to why theoretical rates of 
period change might be different from those measured. If these 
three stars were undergoing mass loss that is about one or two or¬ 
ders of magnitude greater than predicted then the models would 
have rates of period change similar to that observed. However, 
there is no reason why mass loss would be enhanced by such an 
amount for these stars but not the other five, including CMa 
and 6 Cet. 


The goal of our work was to revisit the results of Eggleton & 
Percy! • 1973)1 comparing stellar evolution models and predicted 


rates of period change with observations. Eggleton & Perc^ 
( |1973| l found that stellar evolution models fit observed rates of 
period change for a sample of p Cephei stars well and could ex¬ 
plain the large range of observed rates in terms of both hydrogen- 
core and hydrogen-shell burning stars. While this evolutionary 
scenario explains the diversity of period change measurements in 
P Cephei variable stars, stellar evolution models have advanced 
significantly in the past forty years. 

Since the work of Eggleton & Percy ( 1973| l, stellar opacities 
have changed and new physical processes have been added to 
stellar evolution models. We computed grids of models at solar 
composition, for masses from M - 1 - 20 Mq. The models in¬ 
cluded various amounts of convective core overshooting along 
with different amounts of initial rotation. Eor every model, we 
compute relative rates of period change when the stars evolve 
along the p Cephei instability strip. Both rotation and overshoot¬ 
ing extend the main sequence lifetimes of stars, hence act to de¬ 
crease the relative rates of period change. As such, comparing 
predicted rates of period change to observations do not allow us 
to measure the rotation history of a star or the amounts of over¬ 
shooting because it is not possible to disentangle the stellar mass 
from the amount of mixing in the stellar evolution models. 

We compared our predicted rates to period change measure¬ 
ments for a sample of eight stars ( |Jerzykiewicz|1999 1, but found 
mixed results. Eor three of the stars: y Peg, V2052 Oph and the 
prototype p Cep, our models could not be applied to constrain 
any fundamental stellar parameters. The rates of period change 
for these stars are all consistent with being constant and there 


We computed rates of period change by differentiating the 
period-mean density relation. If p Cephei variable star pulsation 
followed a different relation then the predicted rates of period 
change could be significantly different, but again, if we use a 
different relation then the agreement found for CMa and S Cet 
would be lost. 

A third option is that the measured rates of period change 
for these stars do not represent only secular period change, but 
include the effects of possible period drift where the pulsation of 
the star changes abruptly over short timescales due to some non¬ 
evolutionary cause. Odell| (2012l analyzed timing measurements 
for BW Vul and argued that the pulsation period is mostly con¬ 
stant with respect to time, but underwent random fluctuations 
during various times in the past. These changes are not secu¬ 
lar, i.e., not a function of evolution, but this cause is unclear. A 
similar period drift has been suspected for the classical Cepheid 
Polaris ( [Turner et al.|2005[|Neilson et al.||2012a| ). If this period 
drift is occurring in p Cephei variable stars then it could explain 
why there are such large rates of period change for only three 
stars. It could also explain the two p Cephei variable stars for 
which the rates of period change are well-fit, the stars have not 
undergone period drifts over the time scales of the observations. 

Our results suggest the evolution of p Cephei variable stars 
and their pulsation properties are more complex than expected, 
raising questions about their rates of period change. However, 
our results are based on measurements for only eight stars, ne¬ 
cessitating period change measurements for more p Cephei vari¬ 
able stars along with continued observations of the eight stars 
discussed in Jerzykiewicz ( 1999} . Similarly, our models assume 
that the period-mean density relation is constant and is a valid 
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measure of secular period change that should be verified by di- Turner, D. G., Savoy, J., Derrah, J., Abdel-Sabour Abdel-Latif, M., & Berdnikov, 
rect modelling of pulsation and evolution. L- N. 2005, PASP, 117,207 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295 
Yoon, S.-C. & Langer, N. 2005, A&A, 443, 643 
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